ON THE DUAL VARIABLE OF THE CAUCHY STRESS TENSOR 
IN ISOTROPIC FINITE HYPERELASTICITY 



CLAUDE VALLEE, DANIELLE FORTUNE, AND CAMELIA LERINTIU 



Abstract. Elastic materials are governed by a constitutive law linking the 
second Piola-KirchhofT stress tensor S and the right Cauchy-Green strain ten- 
sor C = F F. Isotropic elastic materials are the special ones for which the 
Cauchy stress tensor a depends solely of the left Cauchy-Green strain tensor 
B = FF T . In this paper we revisit the following property of isotropic hy- 
perelastic materials: if the constitutive law linking £ and C derives from a 
potential a, then a and ln£? are linked by a constitutive law deriving from 
the potential a o exp. We give a new and concise proof which is based on an 
explicit formula expressing the derivative of the exponential of a tensor. 



1. Introduction 

According to the mass conservation principle, the mass density per unit volume 
p and its initial value po are in the ratio 

^ =dctF= (detC)i = (detB)i. 
P 

The relation 

(det F)a = FT,F T 

between the Cauchy stress tensor a and the second Piola-KirchhotT stress tensor £ 
can be rewritten 

- = F — F . 

P Po 

Let us agree to formulate the elastic materials constitutive laws as: 

- = h{C). 
Pa 

The polar decomposition F = RU of the deformation gradient [3] implies: 

B = FF T = RU 2 R T = RCR T or C = R T BR. 
This allows to translate the relation between S and C by a law satisfied by a: 

- = RUh(R T BR)UR T = (RUR T )Rh(R T BR)R T (RUR T ) 
P 

where we have enlightened the tensor RU R T which is nothing else than the square 
root B2 of the positive definite symmetric tensor B. A priori, for elastics materials, 

the tensor — is a function of B and R: 

p 

- = B^ Rh(R T BR)R T B^ . 
P 

It will depend solely of B in a single case: when the tensor Rh(R T BR)R T does 
not depend on the rotation R. The rotations forming a group, the only possible 

l 
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tcnsorial functions h are those satisfying the relations of isotropy with respect to 
B: 

Rh(R T BR)R T = h{B) or R T h(B)R = h(R T BR). 
Because of the relation C — R T BR, the isotropy of the function h can alternatively 
be expressed with respect to C: 

Rh(C)R T = h(RCR T ) or R T h(RCR T )R — h(C). 

To summarize: if the law = h(C) is isotropic, then ^ depends only of B, and it 
is the sole case ; furthermore, under this isotropy condition 

- = B^h(B)B^. 
P 

In this paper, we revisit the property of isotropic hyperelastic materials for which 
the existence of a potential expressing the constitutive law between and C implies 
the existence of a potential linking ^ and ln£>. 

2. Isotropy of the constitutive law linking - and B 

p 

Let f2 be a rotation, if we change B into il T BFl, then B^ is changed into Q T B^fl 
and j is changed in: 

n T B^Qh(n T Bn)n T B^n = n T B^h(B)B^n = n T -n. 

p 

The isotropy of the constitutive law linking and C is thus transferred to the 
constitutive law linking - and B. 

3. COAXIALITY OF B AND h(B) 

Theorem 3.1. Because h is isotropic, the symmetric tensors B and h(B) are 
coaxial (i.e. they have the same eigenvectors). 

Proof. Let n be an eigenvector of B chosen unitary, and let us consider the rotation 
of angle tt around n: 

S = (cos tt)I + (1 — cos ir)nn T — 2nn T — I 

with / as the identity tensor. Such a symmetry S leaves n unchanged and changes 
any orthogonal vector to n in its opposite. The tensor B being symmetric, its other 
two eigenvectors are orthogonal to n, as a consequence S T BS = B. 
The isotropy condition implies S T h(B)S = h(S T BS) or h(B)S = Sh{B), therefore 
h{B)[Sn] = S[h(B)n] or S[h(B)n] — h(B)n. Since the sole vectors unchanged by 
S are the vectors parallel to n, the last equality is possible only when the vector 
h(B)n remains parallel to the vector n, that is to say when n is also an eigenvector 
for h{B). □ 

We easily deduce from this coaxiality property the two next corollaries, which 
will reveal important in the following. 

Corollary 3.2. B and h{B) commute. 

Corollary 3.3. For every real number s, h(B) commutes with the power B s of B. 
The choice s = \ allows one to simplify the expression ^ = B? h(B)B^ in 

- = h(B)B. 
P 
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4. HYPERELASTIC MATERIALS 

4.1. Existence of a potential between the second Piola-Kirchhoff stress 
tensor £ and the right Cauchy-Green strain tensor C. Let us consider a 
derivable function a of C, its derivative Da{C) is a linear mapping from the space 
of symmetric tensors to M. Thus, there exists a symmetric tensor denoted ^ such 
that for every variation 8C of C: 

Da(C)5C = tr(— 5C). 
Hyperelastic materials are those for which there exists a function a such that 

In this assumption, we will say that the constitutive law linking the tensors and 
C derives from the potential a. 

4.2. Derivative of the exponential of a matrix. Let us consider a square ma- 
trix A and a real number t, the exponential exp(tA) is the solution of the matricial 
ordinary differential equation 

— exp(L4) = Aexp(^4) 
dt 

which is equal to / at t = 0. Let 8 A be a variation of A, in the varied equation 

^D(exp)(tA)(tSA) = 8A[exp(tA)] + AD(exp)(tA)(t8A) 

let us introduce the square matrix M(t) defined by 

D(exp)(tA)(t8A) = [exp(tA)]M(t) 
The varied equation becomes 

exp(tA)]M(t) + [exp(tA)}^- = 8A[exp(tA)\ + A[exp[tA)]M (t) 

and simplifies itself into the ordinary differential equation 

^ = [exp(-M)]*A[exp(tA)] 

which can be integrated by quadrature. Because M(0) vanishes, we easily deduce 
from it the value of M(l) and thereafter the variation of the exponential of a matrix 
[10]: 

D(exp)(A)(<L4) = [exp(A)] / [exp(-sA)]8A[exp(sA)]ds. 

Jo 

In the special case where A is the logarithm of the positive definite tensor B, this 
formula allows us to predict that for every variation SB of B: 



D(exp)(lnB)SB = B f B- s SBB s ds 
Jo 
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4.3. Existence of a potential between the Cauchy stress tensor and the 
logarithm of the left Cauchy-Green strain tensor. 

Theorem 4.1. // the tensor derives from a potential a of the tensor C , then 
the tensor — derives from the potential a o exp of the tensor \nB. 

Proof. By deriving the compound function a o exp, we find successively: 
D(a oexp)(ki B)5B = D a (B) (D (exp) (In B) SB) 

o f 1 f 1 

= tr(-^[L>(exp)(m £)££]) = tr(h{B)B \ B- s 8BB s ds) = / tr[h{B)BB- s SBB s )ds. 

To simplify the last integral, it is necessary to pay attention on the switchings 
because the matrix SB does not commute with the others. However, under the 
trace, we can make cross at the beginning the last term of the product of 5 matrices. 
Then from Corollary 2, we can switch this term B s with h(B) and afterwards with 
B, it ends up just before B~ s . The product of the two matrices B s and B~ s reduces 
to the identity tensor /, and the integral simplifies itself into 

ti(h(B)BSB) = tx(-5B). 

P 

The final value of the integral allows to conclude to the constitutive law: 

a d(a o exp) 
~p = a(lnB) 

□ 

5. Conclusion 

Without resorting to the Taylor expansion of the logarithm 11 or of the ex- 
ponential [8] of a symmetric tensor, nor to its spectral decomposition [5], we have 
given an intrinsic proof of the existence of the potential aoexp between — and In B. 
Numerous isotropic hyperelastic constitutive laws expressing directly a in term of 
\nB have been proposed ([2], [B], [7], [S], [H]) and numerically implemented 0]. 
When the potential a o exp is convex, the consideration of its Legendre-Fenchcl- 
Moreau transform is a tool to perform the inversion of the constitutive law ([T], 
[33], [I!])) ie to express the Hencky logarithmic strain tensor lnS in term of the 
Cauchy stress tensor u. 
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